Polarization-based calculation of the dielectric tensor of polar crystals 
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We present a novel method for the calculation of the static and electronic dielectric tensor of 
polar insulating crystals based on concepts from the modern theory of dielectric polarization. As 
an application, we present the first ab initio calculation of the dielectric constants in the wurtzite 
Ill-V nitrides AIN, GaN, and InN. 



The modern quantum theory of polarization in di- 
electrics has been formulated only in recent years [Q. 
This development has opened a new era in the first prin- 
ciples theory of ferroelectricity and pyroelectricity 
The new theory has been used successfully to calculate, 
in a well defined and computationally efficient way, the 
macroscopic polarization changes induced by perturba- 
tions other than an electric field. Examples of such per- 
turbations are e.g. lattice vibrations ferroelectric dis- 
tortions 1^ , and piezoelectric deformations ||^,^ , whereby 
the quantities being calculated are typically dynamical 
Born charges, spontaneous polarization, and piezoelec- 
tric constants. 

So far, no direct attempt has been made towards the 
goal of determining the dielectric tensor (which of course 
quantifies the response to an external electric field) us- 
ing polarization theory. In this Letter we present a novel 
method for calculating the static dielectric tensor of a 
crystal based on concepts from polarization theory; in 
particular, the method rests entirely on the evaluation of 
the dielectric polarization in zero field via the geometric 
quantum phase approach The method works in any 
polar material, i.e. any material having infrared-active 
zone-center modes. As an application, we provide the 
first (to our knowledge) determination of the dielectric 
constants ey along the (0001) axis for the wurtzite III-V 
nitride compounds AIN, GaN, and InN. 

The current method of choice for dielectric response 
calculations is Density Functional Perturbation Theory 
(DFPT) 1^ , a general and powerful approach to response 
properties. The method presented here, besides its differ- 
ent foundations, is less general but considerably simpler 
to implement than DFPT, and it may become a useful 
alternative. 

The static dielectric tensor - The calculation of the 
dielectric tensor is highly non trivial, because it entails 
the determination of the electronic, as well as vibrational 
and elastic-piezoelectric, responses to an external electro- 
static field. The elements of the dielectric tensor are 
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where E is the screened macroscopic electrostatic field 
and P the macroscopic polarization resulting from the re- 
sponse of the electronic and ionic degrees of freedom. In 
the Born-Oppenheimer approximation the macroscopic 
polarization P in the presence of a generic strain or elec- 
trostatic perturbation can be conveniently expressed as 
P = P° -I- P'^* -I- P^, the sum of the spontaneous po- 
larization P'^ of the equilibrium structure in zero field, 
the polarization P'*** induced by lattice response, and 
the electronic screening polarization P^. In the linear 
regime and using Voigt notation, the two latter compo- 
nents can be expressed in terms of the lattice structure 
distortion and screened electric field as 

^" = E^i^' + -^E^"-l (2) 
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where e°° is the electronic component of the dielectric 
tensor, V is the bulk cell volume, ei the strain field, u| 
the displacement of atom s from its equilibrium position; 
e|[^'' is the clamped-ion component of the piezoelectric 
tensor, and Z^f, the Born effective charge: 
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The spontaneous component of the polarization does not 
depend on strain and fields, and it is non- vanishing also 
in zero field. Using Eq. (||) and (||) we can rewrite Eq. 
(Das 
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We shall discuss first how the second and third term are 
evaluated, and then discuss the electronic dielectric ten- 
sor. The last two terms in Eq. (|^) quantify the con- 
tributions of the macroscopic and microscopic structural 
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degrees of freedom to the total polarization, respectively. 
In the absence of other external perturbations, the strain 
field e/ and the atomic displacement are related to the 
screened field E by the condition of vanishing stress, 
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and vanishing HcUmann-Feynman forces, 
Fi = eY,Z:/E, + vY,Etiei-Y,'S>t/i4 =0. (7) 



js' 



The quantities appearing in Eqs. (^) and ^ are the elastic 
constants 
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the piezoelectric tensor 
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and the harmonic force constants and internal strain pa- 
rameters 
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Combining Eq. (|^) and (0), we obtain u and e as a func- 
tion of the electric field: 

< = E '^^k y^Zlf + ^ E -klKrL^^nA E, ■ (11) 

e^^EVfe'efc,^,- (12) 
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These relations are the first two key ingredients of our 
method. Substituting them into Eq. (|^), we obtain after 
some manipulation the following general expression for 
the static dielectric tensor: 
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As implied by Eq. (^, this definition holds for fixed 
stress, certainly an experimentally relevant situation. 
The first term in the last equation is due to atomic dis- 
placements from the ideal position at fixed lattice param- 
eter; the second is the contribution of piezoelectricity- 
related lattice constant changes; the third is the pure 
electronic dielectric screening for a frozen lattice system. 



A central point of the above analysis is that all the 
ingredients of Eqs. (1) and needed for Eq. (|l|) 

can be calculated from distorted and strained bulk cells 
using ab-initio total-energy and force calculations, sup- 
plemented by calculations of the dielectric polarization in 
zero field using the Berry-phase approach |jl|,^||j|]. The 
only exception is the electronic dielectric constant 
for which |6| an alternate approach has to be devised. 

The electronic dielectric tensor - As we now show, 
£°° can be obtained using the relationship between 
macroscopic polarization in zero field and charge accu- 
mulation at the interfaces of an appropriately built ho- 
mojunction of the material of interest In an insu- 

lating superlattice consisting of periodically alternating 
slabs of equal length, stacked along direction n and made 
of materials 1 and 2, the displacement field orthogonal 
to the interfaces is conserved: Di = Ex + Air Pi {Ei ) = 
E2 + ATrP2{E2) = D^- (We use a scalar notation for the 
components of the vectors along n.) Expanding the po- 
larization to first order in the screened fields in the two 
materials as Pi{E) ~ Pi + XiEi, with Pi the polarization 
in zero field and Xi the susceptibility, one obtains 



e^E2. 



In the absence of zero-field (e.g. spontaneous) polariza- 
tion, the familiar equality ef^Ei = £^i?2 is recovered. 
To proceed, we note that periodic boundary conditions 
imply E = Ei = -E2, and /:^E = Ei-E2 = 2E, so that 
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The charge accumulation per unit area at the interface 
between materials 1 and 2 is Sint — ±AE/Att. Therefore, 
switching to an obvious vector notation. 



Si„t = ±2ri-(P2-Pi)/(er+£2 



(14) 



which connects the macroscopic bulk polarizations Pi. 2 
at zero field with the components £^"2 of the dielectric 
tensors of the interfaced materials along the interface nor- 
mal fi. 

In an undistorted honiojunction, i.e. a superlattice in 
which material 1 is identical to material 2, there is no in- 
terface, no polarization change can occur, and the inter- 
face charge is zero. However, a polarization difference can 
be generated in a controlled manner by inducing a small 
distortion S of one of the atomic sublattices in half of the 
superlattice unit cell. The interface charge Sint accumu- 
lated at the interface between distorted and undistorted 
regions can be easily calculated via macroscopic averages 
The zero field polarizations P2 for the material in 
the undistorted state, and Pi for the material in the same 
strain state as in the superlattice, are evaluated directly 
using the Berry phase technique. From Eq. (^), one 
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then extracts the average electronic dielectric constant 
= = +e^)/2. 
In principle Si, the dielectric constant in the distorted 
state, differs from the actual dielectric constant £2; thus, 
so does e. However, in the limit of zero distortion, e 
equals the component of the dielectric tensor along n: 

= Um e. 

This limit can be evaluated with essentially arbitrary ac- 
curacy by extrapolation or interpolation. The procedure 
just outlined, yielding the electronic dielectric constants, 
is the third key ingredient of the method. 

In summary, in the present approach the static dielec- 
tric constant is calculated via (i) calculation of the elastic 
and force constants, (ii) calculation of the piezoelectric 
tensor and Born charges, (iii) evaluation of the electronic 
dielectric tensor. Task (i) requires standard total energy 
and stress calculations; tasks (ii) and (iii) use the geo- 
metric quantum phase polarization; task (iii) also uses 
relatively small, accurately controllable supercell calcu- 
lations. In the latter, one must take care that (a) the 
slabs be short enough that the constant electric field will 
not cause metallization, and that (b) the slabs be suffi- 
ciently long to recover bulk-like behavior away from the 
interfaces. Both requirements are generally met even by 
materials with small calculated gaps for sufficiently small 
applied strains. 

The quantities needed to evaluate the dielectric ten- 
sor are usually obtained by means of DFPT PjTc|]. The 
novelty of the present method is in the absence of a per- 
turbative approach, and in the determination of the elec- 
tronic screening and piezoelectric properties using their 
connection with the geometric quantum phase. 

Application to III- V nitrides - We now apply the for- 
malism just developed to the calculation of £33 |a, the 
component of the static dielectric tensor along the (0001) 
axis at fixed lattice constant a, for the wurtzite III-V ni- 
trides AIN, GaN, and InN jl^. Besides serving as a test 
of our theory, this calculation provides to our knowledge 
the first ab initio theoretical prediction of the dielectric 
constant for these materials. 

In the present case, only a few independent elements 
of the tensors described above are needed, namely those 
containing derivatives of the total energy and polariza- 
tion with respect to the lattice constant c and the in- 
ternal structure parameter u. The piezoelectric constant 
and the Born effective charge involved are 
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The force constant $33 (whence $33^ = 1/<1'33 is ob- 
tained) and the internal strain parameter S33 are cal- 
culated as derivatives of the Hellmann-Feynman force ^3 
with respect to an atomic displacement from equilibrium. 



and to a homogeneous strain of the lattice structure, re- 
spectively: 
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The relevant inverse elastic constant is A 
where 
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All calculations are done using density functional the- 
ory in the local density approximation (LDA) to de- 
scribe the exchange and correlation energy, and ultrasoft 
pseudopotentials for the electron-ion interaction. A 
plane-wave basis cut off at 25 Ry and 12-point Chadi- 
Cohen mesh are found to give fully converged values 
for the bulk properties. Given their known importance 
(Q , the semicore d states of Ga and In are included in the 
valence. The piezoelectric tensor and Born charges have 
been calculated [Q via the Berry phase technique 
using a 16-point Monkhorst-Pack [ p^ k-point mesh in 
the a-plane direction and 10 point uniform mesh in the c 
direction, testing convergence up to 360 total points. For 
the £°° supercell calculation, we have employed (0001)- 
oriented superlattices of typically 4 formula units (16 
atoms), and typical cation sublattice displacements along 
(0001) of 1-2 % of the bond length for AIN and GaN. 
Smaller displacements (~ 0.3 %) were used for InN. 25 
Ry cutoff and 12 k-points in the irreducible super lattice 
Brillouin zone guarantee convergence. No ionic relax- 
ation is allowed, so that the response is purely electronic. 

We report in Table || the calculated dielectric constants 
and their various components as given by Eq. (|l^) for 
AIN, GaN, and InN, along with the available experimen- 
tal data ll6| - [l^ . The electronic dielectric constants are 
very close to the experimental values for both AIN and 
GaN, slightly larger for InN. The calculated static con- 
stant agrees well with experiment for GaN, the only one 
for which it is available experimentally. We also list the 
values of the various constants contributing to £° in Table 
||. For AIN and GaN, the structural constants are simi- 
lar, while the piezoelectric coefficient 633 is much higher 
in AIN. InN behaves somewhat differently, as expected 
from previous experience on other In compounds |^]. 

Discussion - In homopolar semiconductors, the ex- 
ternal field does not cause distortions of the crystal lat- 
tice, so that the static dielectric constant coincides with 
the electronic one; our method as outlined above is not 
applicable in this case, since no zero-field polarization 
can exist in these materials. It does apply, however, to all 
heteropolar materials, in which a polarization can always 
be induced by appropriate atomic displacements. Apart 
from the electronic response, a lower crystal symmetry 
such as in wurtzites enables the action of screening mech- 
anisms related to lattice distortions. As apparent from 
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of Eq. (p^, the difference among e° and e°° is due to 
the polarization induced by optically-active lattice vibra- 
tions (as quantified by the dynamical charges), and to the 
piezoelectric response, if any, along the appropriate axis 
(as measured by the piezoelectric constants). In a previ- 
ous work Q we have shown that AIN, GaN, and InN have 
large effective charges, and the highest piezoelectric co- 
efficients amongst all tetrahedrally bonded semiconduc- 
tors. Indeed (see Table |), the difference between static 
and high-frequency dielectric constants is large: in III- 
V nitrides the phonon-related term is comparable to the 
electronic one, and the piezoelectric component is about 
10% of each of the other two. An important point to be 
noted is that the piezoelectric contribution cannot be ne- 
glected in an accurate calculation. This will be even more 
important in materials with large piezoelectric constants, 
such as ferroelectric perovskites. 

In conclusion, we have presented a novel procedure for 
the calculation of the dielectric tensor based on the geo- 
metric quantum phase polarization theory. The method 
only uses bulk calculation, with the exception of a small 
supercell calculation, needed in the determination of the 
electronic dielectric constant. As an application, we have 
provided the first ab initio prediction of the dielectric 
constants of wurtzite AIN, GaN, and InN. 
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5.06 


0.64 
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GaN 


4.44 


0.15 


5.69 


(5.70*) 


10.28 (10.4^=) 


InN 


5.51 


0.61 


8.49 


(8.40*) 


14.61 



'^From Ref. [0, ^Ref. and '^Ref. (T^ 



(k, k') = Im[ln(detS<^'(k, k'))] 
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TABLE II. Constants needed in the evaluation of e° in 



AIN, GaN, and InN (see text). 





633 


7* 
■^33 


$33 


H33 


A33 


units 


C/m^ 




N/m 


10^" N/m^ 


10''N/m^ 


AIN 


+1.462 


-2.70 


+204 


+0.45 


+3.81 


GaN 


+0.727 


2.72 


+209 


+0.38 


+3.80 


InN 


+1.092 


-3.02 


+155 


+0.32 


+2.22 



